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We propose a new architecture for the measurement-based quantum computation model. The new 
design relies on small composite light-atom primary clusters. These are then assembled into cluster 
arrays using ancillary light modes and the actual computation is run on such a cellular cluster. We 
show how to create the primary clusters, which are Gaussian cluster states composed of both light 
and atomic modes. These are entangled via QND interactions and beamsplitters and the scheme is 
well described within the continuous- variable covariance matrix formalism. 



I. INTRODUCTION 

The one-way model or Measurement Based Quantum 
Computation (MBQC) has emerged as a conceptually in- 
teresting and potentially practical alternative to the stan- 
dard model of quantum computation [T]. MBQC replaces 
the need for coherent unitary control by a sequence of 
adaptive local measurements performed on a highly en- 
tangled resource state known as a cluster state [3 . This 
resource acts as a universal substrate with quantum in- 
formation encoded virtually within it. Though originally 
based on qubits, the cluster model has been generalized 
to higher dimensional discrete-variable systems (qudits) 
|4] as well as to continuous quantum variables [5]. 

In its continuous-variable (CV) incarnation, the re- 
source state is a multimode squeezed Gaussian state. In 
the optical setting [6 , 7 , homodyne detection and photon 
counting plus classical feedforward, suffice to implement 
universal QC over CV. It has also been shown that homo- 
dyne detection alone is sufficient to implement all mul- 
timode Gaussian operations, given a cluster state with a 
sufficiently connected graph [8]. 

Many proposals for the construction of optical CV clus- 
ter states have been put forward, including the linear 
optical construction [6] , the single optical parametric os- 
cillator (OPO) method [S'.TO] and single quantum nonde- 
molition (QND)-gate schemes [TT1[T2]. These procedures 
are deterministic and have the advantage over their dis- 
crete counterparts p^SHTS] . which rely on nondetermin- 
istic interactions and postselection. However, as with 
any CV system, these protocols suffer the usual problems 
such as the lack of infinitely squeezed resources leading 
to finite squeezing errors. 

Here we will use the canonical generation method 
where cluster states are created from single-mode squeez- 
ers and controlled-Z (Cz) gates [5]. The controlled-Z is 
an example of a QND interaction and, for optical modes, 
can be implemented using beamsplitters and inline or 
offline squeezers pIBHIS] . This method, while experimen- 
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tally challenging is achievable with current technology 

m- 

CV cluster states can be built, not only from optical 
modes but also from ensembles of polarized atoms where 
each ensemble is a different CV mode [20^-^3^. The en- 
sembles are entangled by performing interactions with 
off-resonant linearly polarized light, which, on the classi- 
cal level, performs a Faraday rotation. The rotated light 
then serves as a carrier to encode information in each 
of the ensembles and homodyne measurements are per- 
formed to complete the protocol and project the state of 
the ensembles into the desired entangled state. 

However, it has been shown that Gaussian cluster 
states suffer from a fundamental problem. Namely, no 
matter what Gaussian local measurements are performed 
on systems distributed on a general graph, transport and 
processing of quantum information is not possible beyond 
a certain influence region, except for exponentially sup- 
pressed corrections [24 . In other words, as the size of 
the cluster increases there is an exponential decay in the 
Gaussian localizable entanglement over the entire clus- 
ter which leads to the corollary that cluster states with 
only Gaussian operations cannot serve as perfect quan- 
tum wires. Hence large cluster states with only Gaussian 
resources available will always suffer from large errors. 

Here we propose an architecture of a one-way com- 
puter that removes the need for large clusters. In our 
scheme, small, single gate-size cluster states composed of 
atomic ensembles are created and the computations are 
carried out over an array of such clusters. Due to their 
proven suitability as information carriers, communication 
between the clusters is accomplished using optical modes. 
The atomic ensembles we use have the advantage that 
they can serve as a short term quantum memory as well 
as a processing device due to their relatively long lived 
storage times. In order to get information from one pro- 
cessor to the next, we entangle a polarized light pulse to 
the ensemble. This light mode in effect becomes part of 
the processor cluster state. When suitable measurements 
are performed on the atomic modes, the information is 
passed onto the light. This light mode can then be cou- 
pled to the next processor and the procedure is repeated 
until the computation is complete. In order to realize this 
scheme, we require a protocol to include a light mode into 



2 



an atomic cluster state. We call such a state a composite 
cluster state that is composed of m atomic modes and 
n light modes. Here we will explicitly demonstrate the 
protocol for the generation of such states. 

This paper is organized as follow. In section II we will 
review the types of interaction we will be using and show 
how the two-mode composite cluster state is constructed. 
In section III, we discuss the protocol in the covariance 
matrix formalism. In section IV we discuss the entan- 
glement properties of the state. In section V we general- 
ize our discussion to multipartite entanglement and show 
how to create composite cluster states of arbitrary size. 
In section VI we discuss an alternative computational ar- 
chitecture for CV MBQC, making use of the multi-mode 
composite cluster states. We conclude in section VII. 



II. QUANTUM NON-DEMOLITION 
INTERACTIONS 



The physical systems of interest to us are atomic en- 
sembles and light beams. The light beams are used both 
as part of the required cluster state and information car- 
riers between nodes of the cluster. 

Each atomic ensemble contains a large number, Nat, 
of non- interacting alkali atoms with individual total an- 
gular momentum F [25 . The ensemble is described by 
its collective angular momentum J = {Jx-, Jy, Jz), where 



Jk = ^Fk^i, k = x,y,z. 



(1) 



All atoms are assumed to be polarized along the x- 
direction, which corresponds to preparing them in a par- 
ticular hyperfine state, |F, mi?). Then, fluctuations in 
the Jx component of collective spin are kept extremely 
small relative to the strong coherent excitation and we 
can just treat Jx as a classical quantity, Jx ^ (Jx) = 
hJx = hNatomF. The orthogonal components of spin 
are unaffected by the x-polarization and quantum fluc- 
tuations around their mean value remain relatively large. 
By taking an appropriate normalization, the orthogonal 
components fulfill the canonical commutation relations, 
[Jy/^/hJx^Jz/^/hJx] = ih. In this canonical form, we 
can identify these variables as the "position" and "mo- 
mentum" of the system, defined in the following way. 



XA 



Jy 



PA 



Jz 



(2) 



Light modes are used as a node of the cluster and for 
generating interactions. Each light mode is taken to be 
out of resonance from any relevant atomic transition and 
linearly polarized along the x- direction. A useful de- 
scription of light is given through the Stokes operators. 



s = {sx-,Sy^Sz) of light polarization given by 



(3) 



The individual components correspond to the difference 
between the number of photons (per unit time) with x 
and y polarization, ±7r/4 linear polarizations and the two 
circular polarizations, respectively. These allow for a mi- 
croscopic description of the interaction with atoms, how- 
ever only the macroscopic observables Sk = Jq Sk{t)dt^ 
where T is the duration of the light pulse will be rel- 
evant. Similar to the atomic case, the linear polar- 
ization along X allows us to make the approximation 
Sx ^ (Sx) = Nphh/2. The orthogonal components Sy 
and Sz are rescaled to fulfill the commutation relation 
[Sy I V hSx ,Sz/y/ hSx] = ih. Once again we make a con- 
nection with the canonical position and momenta. 



XL 



PL 



(4) 



Note, in the language of the canonical variables, the spin 
and polarization degrees of freedom are treated on an 
equal footing. In the following analysis, we will only use 
these canonical variables to refer to the atomic ensem- 
bles and light modes We denote by Ai{Li) the i^^ atomic 
ensemble (light mode). Let us assume that a light beam 
propagates in the YZ plane and passes through a single 
ensemble at an angle a with respect to the z-direction. 
Then the atom-light interaction can be approximated by 
the effective QND Hamiltonian |26j, 



Heff{a) = — PL (pa cos a + sin a). 



(5) 



where is the coupling constant. The evolution associ- 
ated with this interaction can be evaluated by applying 
the Heisenberg equation for the atoms and the Maxwell- 
Bloch equation (neglecting retardation) for the light. The 
variables characterizing the composite system transform 
according to the following rules: 



Pa — Pa 



pT=Pl- 



- tzpL COS a, 
hzp'^L COS a, 

- hz{p^A COS a 



- x^2 sin a). 



(6) 



These are straightforwardly generalized to the case in 
which a single light beam propagates through many 
atomic ensembles, impinging on the i^^ sample at an an- 
gle ai. 

We also require an interaction between the cluster light 
mode and the interaction light pulses. These interac- 
tions can be implemented by a beamsplitter and squeez- 
ers p!6|, yielding an interaction Hamiltonian of the form 
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H = xlXi^ where xi is the position quadrature of the in- 
teraction hght pulse. Under the influence of this Hamil- 
tonian the variables transform according to: 



^out 



Pi 



Pi 



(7) 



Note that Pi (xi) are the squeezed (antisqueezed) modes 

such that, pL/i = ^~^PL]i ^L/i 

uum modes labelled by the superscript (0). 



e^^x^j]. with vac- 

L/i 



A. Composite cluster protocol 



Hamiltonian H = KXApii- The interaction pulse picks 
up an atomic quadrature term in its position variable, 
^ii ~ ^ • Then the interaction pulse is com- 
bined with the cluster light mode, L, on a beamsplit- 
ter with effective interaction, H = xlXi-^. This modi- 
fies the quadratures of the light mode by rotating to 



p'l = Wl 



KjX , 



Our second interaction pulse i2 




The stabilizer formalism gives us an efficient way to 
represent continuous variable cluster states. A state \(\)) 
is stabilized by an operator K if it is an eigenstate of 
K with unit eigenvalue. That is i^r|^) = |^). If such a 
set of operators exist for a given state, we call that state 
a stabilizer state and we may use the generators of its 
stabilizer group to uniquely specify it. It is well known 
that the stabilizers for continuous variable cluster states 



are 



K,{S)=XAS) n 
ieAr(i) 



l,...,n 



(8) 



for all s G [R, where is the set of vertices that neigh- 
bour vertex i, i.e. = {j\{vi^Vj) G E}. This group 
is described by its Lie algebra, the space of operators H 
such that H\(j)) =0. We refer to any element of this alge- 
bra as a nullifier of |^). Being hermitian, every nullifier 
is an observable [27] . Any ideal cluster state has nullifier 
representation 



Hi Pi 



jeNii) 



i = 1, n. 



(9) 



where Hi ^ for an ideal cluster. 

For the case of the two mode composite cluster com- 
posed from an atomic ensemble and optical mode, these 
nullifiers reduce to just two conditions on the quantum 
variables. 



Pa-xl^ 0, Pl - xa 







(10) 



where the subscript A refers to the atomic ensemble and 
L is the light mode. Using the interaction pulses we can 
assemble these quadrature combinations in the follow- 
ing way. Unlike standard two-mode entangled states, the 
two-mode cluster state mixes the position and momen- 
tum quadratures of different nodes. 

We pass our first interaction pulse, ii though the 
atomic ensemble (see Figjlja)). By careful choice of the 
angle at which the light impinges on the atomic sam- 
ple we can couple the light and the atoms via the QND 



FIG. 1: Schematic description of the creation of a two mode 
composite cluster state. Light modes can be combined with 
atomic ensembles to form the cluster states by enacting a 
QND interaction with squeezers and beamsplitters, (a) The 
first interaction encodes atomic quadratures onto the light 
mode L. (b) The second interaction encodes light quadratures 
onto the atomic ensemble A. These interactions are required 
to fulfil the nullifier condition of Eq. ([22| / 



is first passed through the beamsplitter with light mode 
L (Figjljb)). The quadrature action on 12 is a rotation 
of Pi2, P'i^ = Pl'^ — XL with the associated backaction 
on the light mode. Interaction of the pulse with the 
atomic ensemble via H = nxApi2 results in the output 
pout _ ^pin _ ^^^j _ /^(p^^ — J. Note that the momenta 
of the light and atomic modes now have the correct form 
as given by (22) for the quadrature combinations for a 



two mode cluster state plus some extra noise terms that 
are a consequence of the backaction of the QND interac- 
tions. To complete the protocol, the outgoing interaction 
pulses undergo homodyne measurements which project 
the composite system into a cluster state. As we shall see 
in section VI, such composite states provide the means 
for the atomic clusters to communicate with each other 
so they can act as one computational array. 

Due to the fact that both light and atoms are highly 
polarized, the initial states of the atomic ensembles as 
well as the light can be treated as Gaussian modes. The 
interaction is a bilinear coupling between Stokes opera- 
tors and collective spin and thus preserves the Gaussian 
character of the initial states. Furthermore, the beam- 
splitter interactions that we use to include light modes in 
the cluster also yield a description in terms of Gaussian 
operations. To tackle these CV interactions we employ 
the covariance matrix (CM) formalism. This allows us 
to examine the entanglement properties of the resulting 
state and estimate bounds on the experimental imple- 
mentation. 
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III. SYMPLECTIC DESCRIPTION 

Gaussian functions are mathematically completely de- 
fined by their first and second moments [28l [29l l3T] . 
Hence it follows that any Gaussian state p is charac- 
terized by the first and second moments of the quadra- 
ture field operators. We denote the first moments by 
R = ((i^i), (^2)7 •••7 (Rn)^ i^n)) and the second moments 
by the covariance matrix (CM) a of elements 



1 



(11) 



First moments can be arbitrarily adjusted by local uni- 
tary operations, that is, displacements in phase space. 
These can be performed by applying single-mode Weyl 
operators to re-center the reduced Gaussian correspond- 
ing to each single mode. Such operations leave the struc- 
ture of the Gaussian state and hence any information con- 
tained within it invariant and hence any first moments 
are unimportant to our analysis and from now on, unless 
otherwise stated we will set these to zero. 

The Wigner function of a Gaussian state can be written 
as follows in terms of phase-space quadrature variables 



W{R) = 



(12) 



where R is the real phase-space vector f Ar,PAr). 

A useful observation is that even though the Hilbert 
space in which the state lives is infinite dimensional, a 
complete description of an arbitrary Gaussian state is 
therefore encoded in the 2N x 2 CM a matrix, which we 
will use to denote the second moments of the Gaussian 
state or the state itself. In the formalism of statistical 
mechanics, the CM elements are the two-point truncated 
correlation functions between the 2N canonical continu- 
ous variables. The covariance matrix corresponding to a 
quantum state must fulfill the positivity condition 



i^N > 



where 



N 



N 



0n, n = 



1 
-1 



(13) 



(14) 



In general, for QND and beamsplitter interactions, the 
covariance matrix takes the form 



A C 
B 



(15) 



where the submatrix A corresponds to those modes that 
are nodes of the cluster, B are the light modes responsible 
for mediating the interactions between the nodes within 
the cluster states and C is the correlations between them. 

If a Gaussian state undergos a physical unitary evolu- 
tion that preserves its Gaussianity, then the transforma- 
tion at the level of the CM corresponds to a symplectic 
transformation by a symplectic matrix S 



^out — S O'inS. 



(16) 



We also require the ability to perform homodyne de- 
tections on the outgoing interaction pulses. In the CM 
formalism, assuming an initial displacement of zero, the 
measurement of quadrature xl with outcome zl leaves 
the system in the state described by the CM 



A' = A-C{XBX)-^C^ 
with the displacement 

dA=C{XBX)-\zL,^), 



(17) 



(18) 



where the inverse is understood as the Moore-Penrose 
pseudo-inverse whenever the matrix is not full rank and 
X is a diagonal matrix with the same dimension as B 
with diagonal entries (1, 0, 1, 0, 1, 0). 

We assume the initial state of the composite system is 
given by the covariance matrix for the atoms, light and 
interaction pulses, ain = IL2 ® ^2^ ® where the 
2x2 identity matrices stand for single modes. 

As described in our protocol above, we begin by pass- 
ing ii through the atomic ensemble and then combine 
it with the cluster light mode via a beamsplitter. The 
symplectic matrix for this operation is given by 



Sin 



/ 1 

1 

n 1 




0/^0 



Vo 



\ 
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1 
1/ 



(19) 



The second round of interactions with ^2, is given by the 
symplectic matrix 



Si 



int2 



/ 1 

1 

1 







K 

Vo 



-A^ 

1 





1 
1 
1 




-1 



(20) 



Then the resulting state is described by the transforma- 
tion aout = Sj^t^Sj^^^(JinSint^Sint2 (sec appendix for ex- 
plicit form). However, entanglement between the light 
mode and atomic ensemble is not produced until the in- 
teraction pulses are measured. We perform a homodyne 
measurement on the outgoing pulses ii and 12 in the x- 
basis with results zi and Z2 to project into the final state 
described by the CM 



fin = CF^ — CF^ {Xg^ X) 



(21) 



where cr^ is the upper left 4x4 matrix corresponding to 
the atomic and light modes, cr^ is the lower right 4 ma- 
trix representing the interaction pulses and cr^ and cr^ 
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are the off diagonal matrices with entries corresponding 
to the correlations between the light mode, the atomic 
ensemble and the interaction pulses. The final covari- 
ance matrix for the two-mode composite state is given 
in Appendix. A. The final state is independent of the 
measurement outcomes, but they are present in the dis- 
placement vector. 



IV. VERIFICATION OF ENTANGLEMENT 

On completion of the protocol, an analysis if the cor- 
relations induced between the atomic ensemble and light 
must be performed. Through the CM formalism we have 
complete access to all the information we need to ver- 
ify entanglement. A structural separability test, which 
can only be applied when the full CM is available is the 
positive partial transpose (PPT) test [SOj [32l |33]. In 
phase space, any N mode Gaussian state can be trans- 
formed by symplectic operations in its Williamson di- 
agonal form u [34 , such that a = S^vS^ with v = 
diag{z/i, i^i, i^n, i^Ar}. The set {ui} of all positive- de- 
fined eigenvalues of \iVt(j\ constitutes the symplectic spec- 
trum of cr, the elements of which are the symplectic eigen- 
values which must fulfill the conditions Vi > 1 io en- 
sure the positivity of the density matrix associated with 
cr. The symplectic eigenvalues, are determined by 
N symplectic invariants associated with the character- 
istic polynomial of the matrix |zr^cr|. In order to say 
something about entanglement, we recall that the CM's 
PPT is a necessary and sufficient condition of separabil- 
ity for the (M + A/')-mode bisymmetric Gaussian states 
with respect to the M\N bipartition of the modes. This 
also holds for (M + A^)-mode Gaussian states with fully 
degenerate symplectic spectrum. Further, for the case 
of M = 1, PPT is a necessary and sufficient condition 
for separability of all Gaussian states [35l [36] . In phase 
space, partial transposition corresponds to partial time 
reversal of the CM, or simply, a change of sign of the mo- 
menta for chosen modes. If {ui} is the symplectic spec- 
trum of the partially transposed CM <j, then a (1 + A^)- 
mode Gaussian state with CM a is separable if and only 
if > 1 Vi If the partially time reversed CM does not 
fulfill the positivity condition, the corresponding state 
is entangled. Computing symplectic eigenvalues for the 
partially transposed CM <jj^^, we find that the state is 
indeed entangled since the smallest eigenvalue of \iVtdfin\ 
is z> = 0.63 for an interaction strength of z^: = 0.8. 



V. MULTIMODE COMPOSITE CLUSTER 
STATES 

Here we extend our analysis to the multipartite case. 
In general, composite clusters can be composed of m 
atomic modes with n light modes. We shall call these 
(m, n)-composite cluster states. Here we give the proto- 
col for a (4, 1) cluster, which will become the basic unit 



for the computational scheme in the next section. 

Our protocol for the (4, l)-composite cluster proceeds 
as follows. We construct a four- mode square cluster 
state from atomic ensembles labelled Ai (Figj2|, where 
i = 1,...,4 [20 . We include a light mode to form the 
composite state by entangling it with one of the atomic 
modes in the cluster. 

The nullifiers for a four-mode square cluster can be 
written as: 

PAi - - 0, PA2 - ^Ai - xa^ 0, 

pAs - XAi - XA^ 0, pA^ - XA2 - XA3 ^ 0. (22) 

To entangle the ensembles we use light pulses labelled. 




FIG. 2: A four-node square cluster state, composed of four 
atomic ensembles Ai, A4 

ii, ...,^4. We make use of QND Hamiltonians to mediate 
the interaction between the ensembles and light pulses. 
Note that for such an off-resonant atom- light interac- 
tion, the Hamiltonians are given in [? ] and are well 
established experimentally. The Hamiltonians we make 
use of are Hi = nxApi and H2 = nxAXi which leads to 
quadrature transformations 



Pa^ =VA-'^xf, 
/xout _ /^in 

pT=Pi-^i'X- (23) 

The protocol to generate atomic cluster states is depicted 
in Figjs] (a) Spin information is picked up from A2 and 
A3 by pulse ii via the interaction Hi. This pulse encodes 
the information onto Ai through the interaction H2. (b) 
Spin information is picked up from Ai and A4 by pulse 
12 via the interaction Hi. This pulse encodes the infor- 
mation onto A2 through the interaction H2. (c) Spin 
information is picked up from Ai and ^4 by pulse is via 
the interaction Hi. This pulse encodes the information 
onto A3 through the interaction H2. (d) Spin informa- 
tion is picked up from A2 and A3 by pulse u via the 
interaction Hi. This pulse encodes the information onto 
A4 through the interaction H2. When this procedure is 
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FIG. 3: Atomic cluster state generation protocol (HD - Ho- 
modyne detection): (a)Liglit pulse ii interacts with ensembles 
A2 and A3, picks up atomic spin information which is encoded 
onto Ai. (b)Light pulse 12 interacts with ensembles Ai and 
A4, picks up atomic spin information which is encoded onto 
A2. (c) Light pulse is interacts with ensembles Ai and A4, 
picks up atomic spin information which is encoded onto A3. 
(d) Light pulse 24 interacts with ensembles A2 and A3, picks 
up atomic spin information which is encoded onto A4. The 
interaction pulses are homodyne measured once these steps 
have been completed to form the cluster state. 



completed, the ensembles have quadrature combinations: 





= PAr 


— hi 


— K^^Xa^ 




^Ai 


xa 




= PA2 


— Hi^XA^ 


— Hj^XAx 




^A2 


XA 


^out 
PAs 


= PAs 


— hi^XA^ 


— Hj^XAx 






= XA 


PA4 


= PA^ 


K/ X A2 


Kj X Aq 


- f^N4, 


^A4 


= XA 



(24) 



Pii + Pu^ 
are the 



where A^i = Xi^ ^ ^ Pis, N2 = xi^ 

= Xi^ + + pi^ and 7V4 = Xi^ + Pi^ + Pi^ 
backaction terms due to the QND interactions with the 
interaction pulses. Then for = 1 the remaining terms 
are exactly the nullifier relations, (22), for the four- mode 



square cluster state. To complete the protocol, homo- 
dyne measurements are made on the outgoing interac- 
tion pulses which project the ensembles into the required 
state. 

The (4, l)-composite cluster, has nullifier relations 



where = Xi^ -\-Pi2 +^23 +^25 +Pi6 is the new backaction 
term. The light mode quadratures have transformed as 



pT' =Pl-^^xa,-^Nl, 



(27) 



and Nl = Xi, 



. Note that the backaction terms, 





FIG. 4: To add a light mode to the atomic cluster to form 
the (4, l)-composite cluster, we use the protocol of section II. 
Pulses is and ie interact with ensemble A4 and the light mode 
L via a beamsplitter and are subsequently measured. 



Ni^ are composed of quadratures of the interaction light 
modes only. The interaction modes are momentum- 
squeezed and interact weakly with the ensembles so their 
backaction can be neglected (see also the related experi- 
ment ^25^ where this has been verified). We observe that 
we have a complete set of quadrature combinations that 
satisfies the composite cluster nullifier conditions (25). 
Finally, homodyne measuring and ie completes the 
protocol. 



We note that this protocol can be simply extended 
to creating general (m, n)-composite clusters of arbitrary 
shape. Further QND interactions can be used to add 
atomic ensembles and it is always possible to add a light 
mode to an atomic mode through the beamsplitter in- 
teraction. However adding extra nodes always results in 
extra backaction terms, Nj. 



PAi - XA2 - XA3 0, PA2 - xai - xa^ 0, 

PAs - XAi - XA^ 0, PA^ - XA2 - XAs - XL ^ 

PL - XA, ^ 0. (25) 

We entangle a light mode L with ensemble A4 (FigQ. 
Following the protocol given in section II, the light mode 
is coupled to the existing atomic cluster with interaction 
pulses ^5 and zq. The nullifiers of the atomic modes Ai, 
A2 and As are unaffected but the quadrature combina- 
tions for are now 

pT, = PA, - I^^XA2 - I^^XAs - i^'^xl - K.N^, (26) 



VI. A NEW ARCHITECTURE 

Here we seek to address one of the major difficulties 
in the practical implementation of cluster state compu- 
tation: That of decreasing localizable entanglement with 
increasing size of the cluster when only Gaussian opera- 
tions are available [24 . This places a limit on the size of 
useful cluster states, i.e. those that have sufficient entan- 
glement available to perform processes below the compu- 
tational error threshold. Traditionally, the MBQC model 
relies on creating large clusters on which the entire com- 
putation can be performed and read out. Here we pro- 
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pose to build up computational arrays from small build- 
ing blocks or qubricks which are composed of an atomic 
cluster state acting as a quantum processor and a light 
mode that allow for communication to other qubricks 
(Figj5|. This allows us to eliminate redundant nodes 
that only serve to increase the error rate in the system. 
A similar suggestion to increase the efficiency of gener- 
ating qubit cluster states based on ancilla light modes 
was made in [41 . In our scheme, many atomic cluster 



nication between the processors, we entangle an ancillary 
light mode L with the output node in the atomic clus- 
ter to form a qubrick which is a (4, l)-composite cluster. 
The light mode of the qubrick can then be used to con- 
vey the output of the process performed on the atomic 
cluster contained in the brick to the next in the array 
(Figj6|. Given sufficient resources, the state can be pro- 
cessed in parallel using strings of qubricks. The total 
output from each of the strings can then be combined to 
give the final output. In more detail, a typical process 



Input mode 




Composite cluster 
state 



Output mode 

Atomic ensembles ^ Light mode ^ 



FIG. 5: A basic qubrick composed of a four- mode atomic 
cluster state plus a light mode to create a composite cluster. 

states, labelled (Ci, Cat), are created. Each of these 
act as a small quantum processor which performs some 
unitary operation Ui on an input state. Each processor 
Ci, shares no entanglement with any other in the ar- 
ray. We keep the number of nodes in each atomic cluster 
small, say four nodes each (Figj2|, which is minimum suf- 
ficient for a controlled quantum gate. By limiting the size 
of the cluster states, the decay of entanglement within 
each processor is kept to a minimum. To allow commu- 



l^> 




tl t2 ts 



FIG. 6: The qubricks interact by matching up their output 
light modes with the input of another qubrick. A state can be 
loaded onto the first brick at time ti, the state is processed 
through adaptive measurements and the output encodes onto 
the light mode. The qubrick at ti having been measured no 
longer contributes to the state. The light pulse is connected 
with the next in the sequence of qubrick processes t2 and the 
process is repeated. This operation can be run in parallel 
with other sequences of qubricks and the final state is given 
by their combined outputs. 




FIG. 7: Each cluster d is kept separate and only commu- 
nicate through light pulses Li. The state to be computed is 
loaded onto Ci, which undergos a sequence of adaptive mea- 
surements. A composite cluster is formed to transfer the state 
to a light mode which is carried to the next cluster in the pro- 
cess. Each cluster acts on the state by some unitary Ui and 
the final state is given by lJ\ip). 

would proceed as follows. A state IV^ini), is loaded onto 
the first cluster Ci . Adaptive measurements are then ap- 
plied to the atomic ensembles to process the state with 
outcome IV^O = ^^ilV^im)- A light mode Li, is added to 
Ci to form the qubrick and the state is transferred to 
Li by a further measurement (FigjTj). Since the cluster 
and its associated entanglement has been destroyed in 
the measurement process the light pulse is free to carry 
the information to the next qubrick, without inadver- 
tently entangling the atomic clusters belonging to dif- 
ferent bricks. Another measurement transfers the state 
onto the first node of the atomic cluster, C2 and a gate 
operation is performed by a new sequence of measure- 
ments yielding U2\^|J'). A light mode L2, is added to C2 
to form a new qubrick and it carries the processed state 
is to the next cluster. Since each cluster applies some 
unitary transformation Ui to the state, this process can 
be repeated until the desired output if achieved and the 
state is given by U\iIj) = Ui <S) Ui-i (8) ... <S)U2^ Uil^p). 

It is important to note that we do not claim that this 
procedure outperforms the error rate of large cluster ar- 
rays in general. What the qubrick scheme does guaran- 
tee, is that the loss of entanglement and therefore the 
errors accumulated in each quantum gate are constant. 
Then the errors that propagate through the computation 
are just proportional to the number of qubricks used. 
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This is in contrast to the large cluster state in which the 
entanglement decays exponentially with the size of the 
cluster and hence errors accumulate rather quickly as the 
cluster size increases. Note however, that the error rate 
of particular geometries of large scale cluster states ac- 
tually outperform the qubrick scheme i.e., for a 16- node 
square cluster. In this case the entanglement scales as 
Ce~^, where C is some constant. If we assume the error 
rate increases proportional to the degradation of entan- 
glement then the errors scale as ae^ where a is some con- 
stant that depends on the particular system in question. 
In the qubrick setting, to replicate this state exactly, we 
require four qubricks and the cumulative error then scales 
as ae^, which is significantly worse. Here we do not pro- 
pose that the qubricks mimic the large scale states ex- 
actly (since they can never beat the fundamental limit 
anyway). In this scheme, each qubrick and hence each 
quantum gate is initialized individually and only when it 
is required (Fig.([6|). These gates are then directly cou- 
pled to the light pulse carrying the output from previous 
processes. Our scheme removes the need for nodes serv- 
ing as quantum wires to convey information around the 
cluster and eliminates any redundant nodes. In doing so 
we supress the exponential losses due to entanglement 
decay and hence reduce the error rate. 

Furthermore, this method can be simplified to just one 
qubrick, if we loop the output from the brick back to the 
input in a similar manner to that proposed in [11]. In 
this scheme, the output from the first computation again 
gives = Ui\iIj). The atomic cluster is re-generated 
and the light mode is fed back into the input node where 
a different series of adaptive measurements is performed. 
This gives U2W) which is processed again until the de- 
sired unitary is enacted. This scheme eliminates the par- 
allel element but it vastly decreases the number of re- 
sources required while maintaining a constant error rate. 
This time the error depends on the number of time the 
state is reused. This type of scheme could from the basis 



of a proof of principle demonstration of CV cluster state 
computation with atomic ensembles since it is achievable 
with currently available technology. 



VII. CONCLUSION 

We have presented a scheme to produce a two-mode 
cluster state composed of an atomic ensemble and light 
pulse, called a composite cluster, via QND interactions 
and beamsplitters. We find that since both the atomic 
ensembles and light are initially Gaussian states and 
only Gaussian operations are performed, the procedure 
is well described by the covariance matrix formalism. 
From the covariance matrix we have confirmed that the 
states are indeed entangled using the PPT separability 
criterion. We have generalized our protocol, and given 
an explicit construction for multipartite composite clus- 
ter states composed of an arbitrary number of atomic 
and light modes. Using our composite clusters, we have 
proposed a new architecture for the one- way computer 
that reduces losses due to decay of localizable entangle- 
ment by preparing many separate clusters and using the 
light modes as information carriers between them. This 
scheme is experimentally feasible with current technol- 
ogy and has the potential for scalability since many of 
the resources can be re-used as many times as required. 
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Appendix A: Explicit covariance matrices for 
two-mode composite cluster state 

Here we give the explicit CM for the two- 
mode composite cluster state. Computing aout = 
Slnt2^Inti^inSintiSint2 yields the CM: 
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However, this has not produced the required entangle- 
ment between the atomic and light mode. Homodyning 
the interaction pulses produces the required correlations. 
In terms of covariance matrices, this amounts to the op- 
eration 

a^' =cT^- cT^{XcT^X)-'cT^" (A2) 



which traces out the interaction light modes ii and 12 to 
leave us with a 4 x 4 matrix representing only the atomic 
and light modes that are part of the cluster: 
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(A3) 



This final CM represents the state (up to a known dis- 2x2 block represents the atomic mode, the lower right 
placement that depends on the measurement outcome) 2 x 2 is the light mode and the off diagonal blocks are 
of the two-mode composite cluster. Here, the upper left the correlations between them. 



